Implementations of Nonadiabatic Geometric Quantum Computation using NMR 
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Recently, geometric phases, which is fault tolerate to certain errors intrinsically due to its geo- 
metric property, are getting considerable attention in quantum computing theoretically. So far, only 
one experiment about adiabatic geometric gate with NMR through Berry phase has been reported. 
However, there are two drawbacks in it. First, the adiabatic condition of Berry phase makes such 
gate very slowly. Second, the extra operation to eliminate the dynamic phase. As we know, geomet- 
ric phase can exist both adiabatic(Berry phase) and nonadiabatic (Aharonov-Anandan phase). In 
this letter, we reports the first experimental realization of nonadiabatic geometric gate with NMR 
through conditional-AA phase. In our experiment the gates can be made faster and more easily, 
and the two drawbacks mentioned above are removed. 



Quantum computers can perform certain tasks much 
more efficiently than classical Turing Machine 10. It is 
well known that controlled two-qubit gate, combined 
with single qubit operations, is a universal gate for quan- 
tum computation This two-qubit gate preserves the 
target qubit for the controlling qubit in certain state, say, 
and flips the target qubit for the controlling qubit in 
the other state, Originally this has been achieved ex- 
perimentally using dynamic method in different physical 
systems[| g|. 

On the other hand, central to the experimental real- 
ization of quantum computer is the construction of fault- 
tolerant quantum logic gates In the quest for a low 
noise quantum computing device, geometric phases]?], |), 
which is fault tolerate to certain errors intrinsically due 
to its geometric property, are getting considerable atten- 
tion in quantum computing theoretically || 10, O, [lj, 
H, |l|, |l5[ [t6l |l7|, |H|. Geometric logical gates based on 
Berry phascp | restricted in adiabatic evolution have been 
proposed JlO|, DJ, and the first adiabatic geometric gate 
was implemented in NMR|^]. However, the experiments 
rely on the adiabatic operations. This is a bit imprac- 
tical because the experimental result is inexact unless 
the Hamiltonian changes extremely slowly in the pro- 
cess. But in fact, everything has to be completed within 
the decoherence time. Besides the adiabatic condition, 
both of the previous proposals require extra operation 
to eliminate the dynamic phase. This extra operation is 
unwanted for a fault tolerate gate because if we can not 
eliminate the dynamic phase exactly, the fault tolerate 
property is weakened. For these reasons, one is tempted 
to set up a new scheme which does not rely on the adia- 
batic condition and which does not involve any accessory 
operation performed to eliminate dynamic phase in the 
whole process. 

Indeed, geometric phase does exist in a non-adiabatic 
process. It was shown by Aharonov and Anandan that 
the geometric phase for a two-level system is only de- 
pendent on the area enclosed by the loop on the Bloch 
sphere[|). In the non-adiabatic case, the path of the state 
evolution is, in general, different from the path of the pa- 




FIG. 1: The path evolution of qubit a when qubit 6 is in [T) 
state. The first pulse (Fig. 2) transformed the initial Hamilto- 
nian 2-kJIz to H = 2tvJ (7"cos0 — 7° sin 6) which creates 
an evolution path on the geodesic circle ABC. After time 
t = 1/(2 J), the second pulse changes the Hamiltonian to 
the form 77 = 2-n J f— 7™ cos 9 — 7™ sin #) which creates an 
evolution path on the geodesic circle CDA. Again after time 
t = 1/ (2J), the third pulse restored the Hamiltonian to the 
initial form 2nJI^. Therefore, qubit a undergoes a cyclic 
evolution through a slice circuit C with angle 9 in projective 
Hilbert (density operator) space, and then the AA geometrical 
phase is simply (3 (C) = m Q, where m = ±= is the magnetic 
quantum number and fi = 4$ is the solid angle subtended by 
the slice circuit. 



rameters in the Hamiltonian. The external field need 
not always follow the evolution path of the state, in con- 
trast to the adiabatic case. Thus, it is possible to let the 
external field be perpendicular to the evolution path in- 
stantaneously so that there is no dynamic phase involved 
in the whole process. Recently, several schemes were 
proposed for nonadiabatic geometric gates in different 
systems [f[9[ |20| , pl[ |2^ , p3fl , but there have been no cor- 
responding experimental demonstrations. Here we show 
the realization of nonadiabatic two-qubit gate through 
conditional-AA phase of one two-level subsystem (qubit) 
controlled by the state of another qubit. 

Experimentally, the controlled two-qubit gate was im- 
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FIG. 2: Pulse sequence used to demonstrate controlled- A A 
phase of the state of qubit a or controlled two-qubit gate. The 
black boxes are pulses oscillating at frequency u> a = u a — nj, 
the flip angles are —9, 29 — n and it — 9 from left to right which 
can be realized by choosing the different pulse duration t and 
pulse power P. Here 9 was selected from to n by 9 = 
n — {0, 1, •• • , 16}. All pulses that oscillated at frequency 
Lu a — it J are hard pulses, each pulse duration is bus. Delay 
times are r = 1/2 J between pair of pulses. Then the time of 
nonadiabatic controlled two-qubit gate is about 4.8ms. Note 
this gate time do not depend on the value of AA phase. 



plemented by the nuclear spins of the 1 H and 13 C atoms 
in a Carbon-13 labeled chloroform molecule, the sin- 
gle 1 H nucleus was used as target qubit a, while the 
13 C nucleus was used as controlled qubit 6, |f ) (||)) de- 
scribes the spin state aligned with (against) an exter- 
nally applied, strong static magnetic field Bo in the 
z direction. The reduced Hamiltonian for this two- 
spin system is, to an excellent approximation, given 
by H = u a I% + u b I h z + 2-irJI%I%, where the first two 
terms describe the free precession of spin a ( X B) and 
b ( 13 C I ) about B with frequencies w a /27r w 500Mhz and 
<x>b/27r w 125Mhz. I® are the angular momentum 
operator in the z direction for a (b). The third term of the 
Hamiltonian describes a scalar spin-spin coupling of the 
two spins with J — 21A.9Hz. The spin-spin relaxation 
times are 0.3s for carbon and 0.4s for proton, respec- 
tively. Realistically, the state of this system is a ther- 
mally equilibrium one, while pure 00 state must be pre- 
pared for most quantum computing tasks. To solve this 
problem, a variety of techniques exist to extract from this 
thermal state just the signal from the molecules in the 00 
statep5, 26, 27], i.e. all the spins aligned in +z direction. 
Here we adopt the method of "spatial averaging" to cre- 
ate this effective pure 00 state [p5|. The initial state of 
spin a is |+) = 4| (|T) + |J.)), the start point of the path 
evolution, which is prepared by a 90° pulse oscillating at 
frequency u a along y-axes. This resonance pulse rotate 
the pure state of qubit a from +z-axes to +x-axes. 

To produce controlled-AA phase of the state, we 
change the oscillating frequency uj a of qubit a from to a to 
Lj a — irJ. Now the Hamiltonian of qubit a can be writ- 
ten as H r a ot = (u a -J a ± 7T j) II = (ttJ ± 7T J) 7°, which 

describes the qubit rotating around +z-axis in rotational 
frame with angular velocity uj a = uj a — nJ. That is, when 
qubit b is in the state |f), H™* = 2nJI z ; when qubit b 
in ||), H™* — 0. Consequently, different states of qubit b 



correspond to different frequencies of qubit a in rotational 
frame. It is just this property that we will use to realize 
'controlled-AA phase shift'. In Fig.l, we show the cyclic 
evolution of 1 H nucleus (qubit a) on the Bloch sphere 
when the 13 C nucleus (qubit 6) is in the state The 
pulse sequences to realize this nonadiabatic cyclic evolu- 
tion was shown in Fig. 2. As we know, in toggling frame, 
the track of Hamiltonian is defined by H = H (0)+i/T(t), 
where H (0) is reduced Hamiltonian of the qubit a 
{{n J ± 7rJ) J°), and H\(t) is Hamiltonian of RF pulses. 
Therefore, by applying the pulse with angle along x- 
axes(first black box in Fig. 2), the Hamiltonian is changed 
to H — 2ttJ cos 9 — Iy sm9) in toggling frame. Sub- 
ject to this Hamiltonian, qubit a undergoes procession 
from initial state ^= (|T) + ||})- After time r = 1/ (2J), 

the initial state becomes \<p (r)) = - — ( J^ + ' (|T) — II})- 
Seen from the resulting state \(p{r)) a , the direction of 
spin a is along -x-axis, besides a global phase factor. 
Now by another pulse with angle 29 — it the Hamiltonian 
of qubit a has the form H — 2-kJ {—I" cos 9 — Iy sin 9) . 
Again, subject to this Hamiltonian, qubit a undergoes 
procession from state \ f{r)) a . After time t = 1/(2J), 

the state becomes \ip(2r)) a = (|t) + ||)). The di- 

rection of spin a returns to +x-axis and a phase —29 
appears. Finally we restore Hamiltonian of qubit a to 
the initial form 2ttJI z by the last pulse with angle it — 9. 
Note that this pulse has no influence on the final state of 
qubit a, i.e. \ip{2r)) a . 



From the whole process described above, the Hamil- 
tonian of qubit a has experienced a cyclic evolution be- 
cause of the cyclic property of the pulse sequence. Cor- 
respondingly, an evolution path of ABCDA on the Bloch 
sphere is produced for qubit o, that is ^= (||) + |J.)) — ^ 

(IT) - ID) ^ ^ (IT) + IT))- Obviously the 



angle 9 is the geometric parameter describing the be- 
havior of the Hamiltonian evolution. This geometric pa- 
rameter also describes the phase difference between the 
initial and final state. Note that the dynamic phases ap- 
pearing in two steps of procession cancel out each other. 
The geometric AA phase is equal to (3 (C) = —29. So 
far we have considered the case in which the state of 
qubit 6 is |t). If qubit b is in |D, the Hamiltonian of 
qubit a is zero and nothing will happen to it. In other 
words, qubit a will preserve itself in this case. Therefore, 
the time evolution operator of the pulse sequence (Fig. 2) 
has the property U (2r) |±) a |T) b = e±^ c ) |±} a |f> 6 and 
U (2t) |±) q \i) b = |±)JD b , where |±) corresponds to 
point A and C, respectively, in the Bloch sphere. Hence 
we can regard qubit 6 as a controlling qubit and qubit 
a as the target qubit; controlled-AA phase of qubit a is 
produced depended on the state of qubit b. In the basis 
of |D and |D, the unitary operator that describes this 
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FIG. 3: Experimental values for the controlled-AA phase 
j3 (C) as a function of solid angle fi = 4$. Experimental points 
are shown as small squares; theoretical values are shown as 
smooth curves. We can see the experiment result fit the the- 
ory quite well; the remaining errors could be due to phasing 
process which is influenced by machine noise and non-ideal 
baseline. Besides, pulse imperfect and relaxation also have 
contributions. 



circle evolutin is 
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In particular, when \(3 (C)| = tt/2, this gate is just the 
C-not gate. 

In order to measure the overall AA phase /3 (C) of 
qubit a we also apply a 90° pulse, oscillating at fre- 
quency uj a along y-axes, to transform qubit b in a co- 
herent superposition of states \t/j(0)) b = (||) + |J.)) 6 
before the cyclic evolution. Therefore after this 
cyclic evolution the final state is |+)„ (IT) + ||}) 6 — 
75 l+)a (It) + e~ l(i{C) |1)) 6 - Here, the unobservable AA 
phase /3 (C) = 29 — ^ft of qubit a transfers to inner 
phase of qubit b (13C in our experiments) which can be 
observed with NMR. To do so, we use the signal of initial 
-j= (|T) + state as a reference one, compared with 



the signal of the final ^ (|f) + e' 1 ^ ||)) 6 state; both 
of them are in-phase doublet but out of phase by /3(C). 
Experimentally, we apply additional phasing factor /3 (C) 
to obtain absorptive lineshape after Fourier transforma- 
tion. 

All experiments are performed at room temperature 
and pressure on Bruker Avance DMX-500 spectrometer 
in Laboratory of Structure Biology, University of Science 
and Technology of China. The experimental results are 
shown in Fig. 3. 

Therefore, we have observed controlled-AA phase or 
implemented a nonadiabatic two-qubit gate. Note that 
though the observation of AA phase has been done in 
NMR with a three- level system [p4|, it cannot be used to 
implement universal two-qubit gate. 



Our experiment resolves two drawbacks of the adia- 
batic geometric computation, namely the slow evolution 
and the need of refocusing to eliminate the dynamical 
phases. Let us now compare the gate time of this nonadi- 
abatic geometric gate to that of adiabatic geometric gate 
and dynamic gate. Since the gate time is limited directly 
by the strength of coupling constant J of the sample. 
Two experiments we selected to compare used the same 
sample as ours[|[ S|, that is, Carbon-13 labelled chloro- 
form sample. In our experiment, it took about 4.8ms to 
realize the gate, slightly longer than the time used to re- 
alize dynamic two-qubit gate (about 2.4ms) ||, yet much 
shorter than the time it took to realize the adiabatic ge- 
ometric two-qubit gate (about 120ms) |)). As the adia- 
batic geometric gate operates for a significantly longer 
time, it is much more severely affected by decoherence. 
This has serious implications for the physical realization 
of adiabatic geometric quantum computation. On the 
other hand, since the state is always perpendicular to 
the effective magnetic field, there is no dynamical phase 
accumulation during the evolution, hence the resulted 
phase factor after cyclic evolution was pure geometric 
phase. Although this nonadiabatic geometric gate is ex- 
perimentally realized in the NMR system, the basic idea 
is general, and could be applied in other physical systems. 
We believe our experiment has led the idea of geometric 
quantum computation much more practical than before. 
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